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Abstract 

J^ ■ We present the full 0{aew) electroweak radiative corrections to the Higgs-boson pro- 

*^ , duction in association with Z'^-boson pair at an electron-positron linear collider(LC) in 

the standard model. We analyze the dependence of the full one-loop corrections on the 
"^ i Higgs-boson mass rriH and colliding energy ^/s. We find that the corrections significantly 

^^ • suppress the Born cross section, and the 0{aew) electroweak radiative corrections are 

generally between 1.0% and —15% in our chosen parameter space, which should be taken 

into consideration in the future precise experiments. 



o 



X 



PACS: 12.15.Lk, 14.80.Bn, 14.70.Hp,11.80.Fv 



'Supported by National Natural Science Foundation of China. 



I Introduction 

One of the most important missions of the future high energy experiments is to search 
for scalar Higgs-boson, which is beheved to be responsible for the breaking of the elec- 
troweak symmetry and the generation of masses for the fundamental particles in the standard 
model(SM)^ [2]. Until now the Higgs-boson hasn't been observed yet, only LEP II experi- 
ment provides a lower bound of 114.4 GeV[H] and an upper bound of 260 GeVj4j for the mass 
of the SM Higgs-boson at 95% confidence level. People believe that with the help of future 
high energy colliders, such as the CERN large hadron collider (LHC), and the linear colliders, 
TESLA, NLC, JLC and CERN CLIC, the existence of the Higgs-boson would be proved or 
excluded in the experiments. 

As far as we know, the present precise experimental data have shown an excellent agree- 
ment with the predictions of the SM except the Higgs sector j^. These data strongly con- 
strain the couplings of the gauge-boson to fermions {gzff aiid g^yrf,), and gauge bosons 
self-couplings, but say little about the couplings between Higgs-boson and gauge bosons, 
which wouldn't exist if the corresponding scalar field has no vacuum expectation value. In 
order to reconstruct Higgs potential, the precise predictions for Higgs couplings, which include 
Yukawa couplings, the couplings of Higgs to gauge bosons and the Higgs self-couplings, are 
necessary. At an e^e~ linear collider with y/s ~ 300 — 500 GeV, Higgs-boson with an interme- 
diate mass value would be produced mainly via the Higgs strahlung process e^e~ — > Z^H^ , 
and the coupling of Higgs-boson to Z-bosons is probed best in the measurement of the cross 
sections of the Higgs strahlung process e^e~ — > Z^H^ and the WW/ZZ fusion processes 
e^e" -^ H^vv and e^e~ -^ H^e'^e". In Refs.jO], it shows that the coupling gzzH can be 
determined at a few percent level for a 120 GeV Higgs-boson with an integrated luminosity 
of 500 /6~^ from the production cross section through the process e^e~ — > Z^H^ . There is 
another class of processes which is interesting for the studies of Higgs physics at linear collider 
called Higgs-boson production in association with a pair of final particles which can be used 



to test the Yukawa couplings, couplings between Higgs-boson and gauge bosons, and Higgs 
self-couplings. For example, e~^e~ -^ Z^Z^H^ process is not only an important process in 
probing gzzHi but also possible to provide further tests for the quadrilinear couplings (such 
as C-violating HZZZ or H^ZZ), which do not exist at tree- level in the SM, because these 
quadrilinear couplings would induce deviations from the SM predicted observablesT^. We be- 
lieve that once the neutral Higgs-boson is discovered and its mass is determined, the double 
Z'^-bosons production through e~^e~ -^ Z^Z^H^ process may provide the detail information 
of the coupling between Higgs-boson and Z^ gauge bosons, which directly reflects the role 
of the Higgs-boson in electroweak symmetry breaking. Moreover, a theoretical accurate esti- 
mate of this class of processes is essential, since e~'"e^ -^ Z^Z^H^ process could be potential 
backgrounds for possible new physics. 

There have been already some theoretical works in investigating the Yukawa couplings, 
the Higgs self-couplings and Higgs couplings with gauge boson pair in the SM at LC, for 
example, the calculations of the NLO QCD and one-loop electroweak corrections to the 
e"'"e~ —>■ tiH^ process in Refs.jHl El ^3 ^| and 77 —>■ ttH^ process in Ref.^21 in probing 
Yukawa coupling, the one-loop electroweak corrections to the process e^e'' -^ H^H^Z^ in 
Ref. \i'A\ I14j for testing Higgs self-coupling. The Higgs productions in association with vector 
gauge bosons(e+e- -^ H^W+W , e+e" -^ H^Z^Z^ and e+e" -^ H^Z^-f) in the SM for 
testing the couplings between Higgs-boson and gauge bosons were studied at the tree-level 
in Ref.[3|. In this work we calculate the full one-loop electroweak corrections to the process 
e^e~ -^ Z^Z^H^ in the SM. The paper is arranged as follows: In Section II we give the 
analytical calculations of the Born cross section and the full 0{aew) electroweak corrections 
to the e~^e~ — > Z^Z^H^ process. In Section III we present some numerical results, and finally 
a short summary is given. 



II Calculation of e+e 
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The tree-level Feynman diagrams contributing to the process e'^e~ -^ Z^Z^H^ in the frame 
of the SM are depicted in Fig^ Due to the fact that the Yukawa coupling strength between 
Higgs/Goldstone and fermions is proportional to the fermion mass, it is reasonable to neglect 



the contributions of the Feynman diagrams which include H^ 
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Figure 1: The tree-level Feynman diagrams for the process e'^e — > Z^Z^H^ . 

We calculated the Born cross section of the process e^e~ -^ Z^Z^H^ by using 't Hooft- 
Feynman gauge and unitary gauge to check the gauge invariance by adopting FeynArts 3.2 
package fro], and got the coincident numerical results. The electroweak one- loop Feynman di- 
agrams can be classified into self-energy, triangle, box and pentagon diagrams. As a represen- 
tative selection, the pentagon diagrams are depicted in Fig. 2. Their corresponding amplitudes 
may involve five point tensor integrals up to rank 4. In the amplitude calculation of the pro- 
cess e^e~ — > Z^Z^H^ involving one-loop contributions, we create all the tree-level, one-loop 
Feynman diagrams and their relevant amplitudes in the 't Hooft-Feynman gauge by using 
FeynArts 3.2, and the Feynman amplitudes are subsequently reduced by FormCalc 4 |16j . 
Our renormalization procedure is implemented in these packages. The numerical calculation 
of the two-, three- and four-point integral functions are done by using FF package J17j. The 
implementations of the scalar and the tensor five-point integrals are done exactly by using 
the Fortran programs as used in our previous works on e^e~ -^ tiH^ and e^e~ -^ Z^H^H^ 
processes [HIE] with the approach presented in Ref . P^ . 




Figure 2: The pentagon Feynman diagrams for the process e~^e -^ Z^Z^H^ 



The 0{aew) virtual electroweak correction to the e^e -^ Z^Z^H^ process can be ex- 
pressed as: 

O-virtual = CTtreef^virtual = 9 OpT^ / '^^3^^e(A^treoA4|ij.,.^^l) (1) 

spin 

where the first factor 2 comes from the two identical Z'^-bosons in the final state, pi is 
the momentum of the incoming positron in the center of mass system(c.m.s.), d<l>3 is the 
three-body phase space element, and the bar over summation recalls averaging over initial 
spins ^nj- o'tree and A^tree are the cross section and amplitude at the tree-level for pro- 
cess e^e~ -^ Z^Z^H^ , respectively. A^virtuai is the renormalized amplitude from all the 
electroweak one-loop Feynman diagrams and the corresponding counterterms. The related 
renormalized quantities and renormalization constants are defined as in Ref. j22j . and can be 
evaluated by using the corresponding equations shown in this reference. 

The total unrenormalized amplitude corresponding to all the one-loop Feynman diagrams 
contains both ultraviolet (UV) and infrared (IR) divergences. To regularize the UV diver- 
gences in loop integrals, we adopt the dimensional regularization scheme |2U| in which the 
dimensions of spinor and space-time manifolds are extended to D = 4 — 2e. And we adopt the 
on-mass-shell (OMS) scheme |211 122| to renormalize the relevant fields. All the tensor coeffi- 
cients of the one-loop integrals can be calculated by using the reduction formulae presented 
in Refs.j2Sl HI]- We check the UV finiteness of our results of the whole contributions of the 
virtual one-loop diagrams and counterterms both analytically and numerically by regularizing 
the IR divergence with a fictitious photon mass. As we expect, the UV divergence contributed 
by virtual one-loop diagrams can be cancelled by that contributed from the counterterms ex- 
actly. 

The soft IR divergence in the process e"'"e~ — > Z^Z^H^ is originated from virtual pho- 
tonic corrections, which can be exactly cancelled by adding the real photonic bremsstrahlung 



corrections to this process in the soft photon hmit. In the real photon emission process 

e+bi) + e-(p2) ^ Z'^ih) + Z0(A:2) + ^^(feg) + l{k^), (2) 

a real photon radiates from the electron/positron, and can have either soft or collinear nature. 
The collinear singularity is regularized by keeping electron(positron) mass. We use the general 
phase-space-slicing (PSS) method 25 to isolate the soft photon emission singularity in the 
real photon emission process. In the PSS method, the bremsstrahlung phase space is divided 
into singular and non-singular regions, and the cross section of the real photon emission 
process (J2| is decomposed into soft and hard terms 

cereal = O"soft + 0"hard = Crtreci^soit + (^hard)- (3) 

where the 'soft' and 'hard' describe the energy of the radiated photon. The energy E^ of 
the radiated photon in the center of mass system(c.m.s.) frame is considered soft and hard if 
E^ < AE and E^ > AE, respectively. Both cisoft and cJhard depend on the arbitrary soft cutoff 
AE/Eb, where Ef, = y/sjl is the electron beam energy in the c.m.s. frame, but the total cross 
section of the real photon emission process (Treai is cutoff independent. Since the soft cutoff 
AE/Eb is taken to be a small value in our calculations, the terms of order AE/E^ can be 
neglected and the soft contribution can be evaluated by using the soft photon approximation 
analytically PHI ^ On ] 

aew f d?k^ ( pi P2 



2vr^ J\k-,\<AE 2^7 \pi -k^ P2-k^J 
As shown in Eq.Q), the soft contribution has an IR singularity at m^ = 0, which can be 
cancelled exactly with that from the virtual photonic corrections. Therefore, i7virtuai+soft; the 
sum of the 0{a^^) virtual and soft photon emission cross section corrections, is independent 
of the fictitious small photon mass m^. The hard contribution, which is UV and IR finite, 
is computed by using the Monte Carlo technique. Finally, the corrected cross section for the 
e^e^ — > Z^Z^H^ process up to the order of ©(af^) is obtained by summing the 0{al^) Born 



cross section utree, the 0{a^^) virtual cross section cTvirtuah and the 0{a^^) cross section of 
the real photon emission process (HJ, i.e., 

<7total = CTtree + 0"virtual + <^Tea\ = <7tree (1 + <^total) , (5) 

In order to analyze the origins of the one-loop electroweak corrections clearly, we organize 
the full one-loop electroweak corrections to the process e+e^ -^ Z^Z^H^ in two parts, the 
QED correction part and the weak correction part. The QED part comes from the QED 
virtual correction and the real photon emission correction. The QED virtual correction in- 
volves the contributions from the one-loop diagrams with virtual photon exchange in the 
loop. For some counterterms involved in the QED contribution, we only have to take into 
account purely photonic contribution to the wave function renormalization constants of the 
electron/positron. The rest of the total virtual electroweak corrections is called the weak 
correction part. With such definitions of the origins of the radiative correction we can divide 
the full one-loop electroweak corrected total cross section in following form: 

_ I I _ I QED I QED I W 

O'total — O'tree + ^virtual+soft + ^hard — 0"tree + "^virtual+soft ' "^hard "'" ^ 

= atree (l + .^^^'^ + S^) = atree (1 + <5total) , (6) 

where (Tvirtuai+soft is the cross section correction contributed by the virtual electroweak one- 
loop diagrams and the soft photon emission process, i^^rtuai+soft' "^hard ^^*^ "^^ ^^^ ^^^ 
corrections from the QED contributions including only the photonic one-loop diagrams and 
the soft photon emission process, the hard photon emission process and the weak virtual 
contribution, separately, (^virtuai+soft ^^ ^^^ relative correction including the contributions of 
QED one-loop diagrams and the soft photon emission process. 6^^^, 6^ and (5totai are the 
relative corrections contributed by the QED correction part, the weak correction part and 
the total electroweak correction, respectively. 



Ill Numerical Results and Discussions 

In our numerical calculation, we adopt the aew-scheme, and the input parameters are taken 
as follows [T^: 

me = 0.510998902 MeV, m^ = 105.658369 MeV, m^ = 1776.99 MeV, 

m„ = 66 MeV, rUc = 1.2 GeV, mt = 178.1 GeV, 

rud = 66 MeV, m^ = 150 MeV, mj, = 4.3 GeV, 

mw = 80.425 GeV, niz = 91.1876 GeV. aew{0) = 1/137.036 (7) 

There we take the electric charge defined in the Thomson limit ae«,(0) = 1/137.036 and the 
effective values of the light quark masses (m^ and rrid) which can reproduce the hadronic 
contribution to the shift in the fine structure constant aew,(w,^) j27j . 

Besides the parameters mentioned above, we must provide the values of the regulator m^ 
and soft cutoff AE/Eb. As we know, the total cross section should have no relation with 
these two parameters. We checked the photon mass independence of the total cross section, 
and found that in the cases of ttij = 10^^^ GeV and ttij = 1 GeV, the 0{a'^^) cross sections 
(Tvirtuai+soft are invariable within the statistical error 1.0 x 10~^. In the numerical calculation 
it shows that if we take a very small value for the regulator rriy, the Monte Carlo sampling of 
the virtual correction part is very slow and it will take a long time to get requested accuracy. 
With this consideration we take m-y = 1 GeV and AE/Ef, = 10^^, if there is no other 
statement. 

To show the independence of the total correction on the soft cutoff AE/Eb, we present 
the 0{aew) relative correction to the e^e~ -^ Z^Z^H^ process as a function of AE/E^ in 
FigEl with niH = 115 GeV and ^/s = 500 GeV. As shown in this figure, both (5soft+virtuai and 
<Jhard obviously depend on the soft cutoff AE/E^, but the full 0{aew) electroweak relative 
correction Jtotai is independent of the soft cutoff value. 





£..\J - 


1 1 1 1 1 1 1 




1.5- 


Vs=500 GeV 
m =115GeV 




1.0- 


---- 5 " 

~ - - _ _ hard 






~ - - - _ _ 




0.5- 


~ ~ - - _ _ 


r^o 


0.0- 
-0.5- 


total 




...-■■••■■" 




-1.0- 


..--•■' virtual+soft 








-1.5- 


- 




-2.0- 


1 1 1 1 1 1 1 



Log,o(AE/EJ 



Figure 3: The 0{aew) relative correction to the e"'"e — > Z^Z^H^ process as a function of 
the soft cutoff l^E/Eb 



In Fig|l]we present the Born cross section cJtree and the full one-loop electroweak corrected 
cross section fitotai as the functions of the c.m.s. energy ^/s with rriH = 115 GeV, 150 GeV 
and 200 GeV, respectively. We can see that the corrected cross sections are always less than 
the corresponding tree- level cross sections clearly, i.e., the radiative corrections are always 
negative in our chosen parameter spaces. We also find that the curves for both (Jtree and Utotai 
go up rapidly to reach their maximal values with the increment of \/s in the region near the 
threshold, and then go down to approach small values. We can read out from the figure that 
the cross sections Utree and Ctotai reach their maximal values of about 0.555 fh and 0.529 fb 
respectively, in the vicinity of ^/s ~ 500 GeV when mu = 115 GeV. For niH = 150 GeV, 
the maximal values of utrcc and (Ttotai are at the position around ^/s ~ 600 GeV and have the 
values about 0.342 fh and 0.320 fb separately. When mn = 200 GeV, the o"tree and citotai 
can reach about 0.200 fb and 0.184 fb at ^/s ~ 700 GeV, respectively. 

FigEldisplays the full 0{aew) electroweak relative correction for the e^e~ -^ Z^Z^H^ process 
versus the c.m.s. energy ^/s, with uih = 115 GeV, 150 GeV and 200 GeV, respectively. As 

10 



shown in the figure, the full 0{aew) electroweak corrections suppress the Born cross sections 

in the range of 300 GeV < \/s < 2000 GeV. The relative correction can be beyond —30% near 

y/s = 300 GeV with mn = 115 GeV, but that value near the threshold is phenomenologically 

insignificant. When the c.m.s. energy is far beyond the threshold, the relative correction 

becomes insensitive to ^/s. The values of the relative corrections vary in the ranges from 

-0.32% to -2.88%, from -4.28% to -7.14%, and from -7.64% to -10.6% when ^ goes up 

from 800 GeV to 2000 GeV for niu = llSGeV, 150 GeV and 200 GeV, respectively. 
0.6- 
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Figure 4: The Born and one-loop level corrected cross sections for the e^e~ 
Z^Z^H^ process as the functions of the e~^e~ colliding energy ^/s 



We plot the Born cross section cjtree and the electroweak corrected cross section cJtotai as 
the functions of the Higgs-boson mass rriH in Fig 1^1 with different values of y/s. It shows that 
both cJtrcc and atotai decrease with the increment of the Higgs-boson mass, and the less the 
value of ^/s, the more rapidly they drop. We can see in this figure that the cross sections of 
(Jtree and (Ttotai with ^/s = 2000 GeV are insensitive to tjih, and the corresponding curves are 
stable as the increment of mn from 100 GeV to 200 GeV. The figure shows also that each 
curve for the one-loop corrected cross section has two spikes, which just reflect the resonance 
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Figure 5: The 0{aew) relative corrections to the e^e~ — > Z^Z^H^ process as the functions 
of ^/s 

effects of the virtual corrections at the positions of mn = 2mw and mn = '^inz, respectively. 
Since we do not take the complex masses of the W^ and Z^ bosons in the calculation of 
one-loop integrals, the numerical results in the vicinities of rriH = 2mw and rriH = 2mz are 
not reliable. 

Figd shows the relationship between the full 0{aew) relative correction of the e~^e~ — > 
Z^Z^H^ process and the Higgs-boson mass tuh- We find that the curves for different val- 
ues of colliding energy decrease with the increment of Higgs-boson mass mn- The relative 
corrections for y/s = 500 GeV have negative values, and are generally smaller than the cor- 
responding ones for y/s = 800 GeV and ^/s = 2000 GeV. We can see from the figure that on 
each curve there are also two spikes at the positions of mn = 2m\\r and mn = 2mz due to the 
resonance effects. Again the numerical relative corrections in those positions are unreliable 
for the same reason as declared for FigEl For ^/s = 500 GeV, the most favorable colliding 
energy for e^e~ — > Z^Z^H^ process with intermediate Higgs-boson mass, the relative correc- 
tion decreases from —3.24% to —13.8% as niH increases from 100 GeV to 200 GeV . We can 
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also read out from the figure that when Higgs-boson mass mn increases from 100 GeV to 
200 GeV, the relative correction for ^/s = 800 GeV goes down from 0.728% to —7.82%, and 
from -1.29% to -11.1% for ^ = 2000 GeV. 







140 160 180 200 

m^ [GeV] 

Figure 6: The Born cross section and the one-loop level corrected cross section as the functions 
of the Higgs-boson mass rriH 

In Table 1 we list and compare some numerical results of the radiative corrections to 
e~^e~ -^ Z^Z^H^ process contributed by the QED correction part, the weak correction part 
and the total electroweak one-loop correction. In our calculation we set the soft cutoff AE/Eh 
being 10~^. From the table we can see that when ^/s = 500 GeV and m// = 115, 150 GeV, 
the full QED one-loop corrections have negative signs and are much larger than the weak 
corrections, while when y^ = 1000 GeV and mn = 115, 150 GeV, the weak corrections 
become important and their absolute values of the weak correction part are comparable with 
those of the QED correction part. For the case with ^/s = 1000 GeV and rriH = 115 GeV, 
the large cancellation between the corrections from the QED and weak parts makes a relative 
small total one-loop electroweak correction. 
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Figure 7: The full 0{aew) relative corrections to the e'^e~ -^ Z^Z^H^ process as the functions 
of the Higgs-boson mass rriH 
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Stotal{%) 


<5«^^(%) 


6^{%) 


500 


115 
150 


-0.2587 
-0.1491 


0.2333 
0.1223 


-0.0002 
-0.0014 


-4.621 

-8.857 


-4.576 
-8.412 


-0.045 
-0.445 


1000 


115 
150 


-0.1538 
-0.1191 


0.1738 
0.1293 


-0.0210 
-0.0207 


-0.320 

-4.317 


6.360 

4.204 


-6.680 
-8.521 



Table 1: The comparison between the numerical results contributed by the QED cor- 
rection part, the weak correction part and the total electroweak 0{aew) correction to 
e'^e" -^ Z^Z^H^ process 
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IV Summary 

In this paper we calculate the full 0{aew) electroweak radiative correction to the e^e~ -^ 
Z^Z^H^ process within the framework of the SM at linear colliders. We analyze the depen- 
dence of the Born cross section, the full 0{aew) electroweak corrected cross section and the 
relative correction on colliding energy ^/s and Higgs-boson mass mn- From the numerical 
results we find that the full 0{aew) electroweak correction significantly suppresses the Born 
cross section. When ^/s = 2000 GeV, both the Born and corrected cross sections are insen- 
sitive to Higgs-boson mass ran- But when ^/s is relatively smaller, e.g., \/s = 500 GeV or 
800 GeV, both the Born and corrected cross sections decrease sharply with the increment 
of Higgs-boson mass mn- In our chosen parameter space, the relative corrections are in the 
value range between 1.0% and —15%. These corrections are so remarkable that we must 
consider them in the precise experimental analysis. 
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Figure Captions 

Figure 1 The tree-level Feynnian diagrams for the process e"'"e~ -^ Z^Z^H^ . 

Figure 2 The pentagon Feynman diagrams for the process e"^e~ -^ Z^Z^H^ . 

Figure 3 The ©(ctew) relative correction to the e~^e~ -^ Z^Z^H^ process as a function 
of the soft cutoff AE/Eb. 

Figure 4 The Born and one-loop level corrected cross sections for the e~^e~ -^ Z^Z^H^ process 
as the functions of the e"'"e~ colliding energy ^/s. 

Figure 5 The 0{aew) relative corrections to the e^e~ -^ Z^Z^H^ process as the functions 
of ^/s. 

Figure 6 The Born cross section and the one-loop level corrected cross section as the 
functions of the Higgs-boson mass mn ■ 

Figure 7 The full 0{aew) relative corrections to the e+e^ -^ Z^Z^H^ process as the 
functions of the Higgs-boson mass niH ■ 
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